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Abstract 

The identification of new materials with superconducting properties is the pursuit 
in the realm of superconductivity research. Here, excitedly, we show that the 
simplest salt daily used can be made a superconductor at normal pressure only by 
adjusting its stoichiometry of Na and Cl as Na3Cl at normal pressure based on 
first-principles calculations. This bulk stable abnormal Na-Cl stoichiometric 
crystal of 3:1, the first ‘magic’ ratio, includes metallic (Na) atoms in the core as 
well as hybridization of ionic and metallic bonding, facilitating the electron- 
phonon-coupling for superconductivity with a critical temperature Tc of 0.13 K. 
The flat bands and van Hove singularities near the Fermi level produce large 
densities of states, similar to H3S and LaHı0, which is beneficial for the emergence 
of superconductivity. The crystal composed of with abnormal Na-Cl magic 
stoichiometry is a precisely tunable, purely sodium and chloride-based, three- 
dimensional bulk superconductor, which is therefore an ideal material for 
designing and understanding abnormal stoichiometric crystals. The methodology 
of constructing this bulk abnormal crystal may be general to almost all elements, 
which could lead to insights into the physics of other conventional superconductors 
and even high-critical-temperature superconductors. 


Superconductors have always obtained great attention!*. However, up to now, it 
remains an open question of whether superconductivity is possessed by most and even 


all elements or some given elements of a specific characteristic. This is the key to 
overcome challenges in the field of superconductivity research to identify new kind of 
materials and search room-temperature superconducting materials at normal pressure. 

Historically, bulk crystals composed of main group and transition metal elements, 
exhibiting unconventional stoichiometry, have predominantly emerged under ultra-high 
pressure conditions*, yielding remarkable properties such as near-room-temperature 
superconductivity>’, transparent inorganic electrides*, and high-energy-density 
characteristics’. However, stabilizing such structures under ambient conditions has long 
posed a formidable challenge!”'’. Based on traditional stoichiometry frameworks to 
find new superconductivity, especially for room-temperature superconducting materials, 
the research has encountered bottlenecks and restricted our thinking. Exotic 
stoichiometric bulk crystals at normal pressure have a great potential to disclose 
exciting properties!*)> , 

Here, by only adjusting the stoichiometry of Na and Cl as the magic ratio Na3Cl, 
the simplest salt daily used can be made a superconductor at normal pressure based on 
first-principles calculations. This cation-rich abnormal bulk Na3Cl crystal includes 
metallic (Na) atoms in the core and hybridization of ionic and metallic bonding, 
facilitating the electron-phonon-coupling for superconductivity, and we show its 
superconductivity with a critical temperature Tc of 0.13 K. The methodology of this 
bulk abnormal crystal construction may be general to almost all elements, illustrating 
that superconductivity maybe possessed by most and even all elements. 


Computational Methods 

We performed the structure searches for the binary (Na3Cl)x system (x = 1—5) using 
the IM*ODE (Inverse Design of Materials by Multi-objective Differential Evolution)’*. 
The calculations of electronic properties were carried out by using PAW potentials in 
the Vienna ab initio simulation package (VASP)!’. The projector augmented waves 
(PAW) approach was used to describe the interaction between ions and electrons!*!”, 
and the exchange-correlation functional was described by adopting the Perdew-Burke- 
Ernzerhof (PBE) generalized gradient approximation (GGA) 7°. We chose the cutoff 
energy of 1000 eV and Monkhorst-Pack k meshes of 2m x 0.03 A’! to ensure the 
enthalpy calculations accurately converge to 1 meV/atom. Lattice-dynamical and 
superconducting properties were estimated using the Allen-Dynes modified McMillan 


2122 with ultrasoft 


equation, performed within the Quantum ESPRESSO package 
pseudopotential adopted. The kinetic energy cutoff was set to 50 Ry. The k-points of 12 


x 12 x | and q-points of 4 x 4 x 1 were employed. 


Results and Discussion 


The ground state structure of Na3Cl found by global optimization without external 


pressure is a six-layer structure with space group P2;/m (number 11), as depicted in 
Figure la and 1b. There are 2 Cl atoms and 6 Na atoms in the unit cell of Na3Cl_ P2)/m, 
and from bottom up, it is composed of Na-Cl layer, 4 pure Na layers, and then Na-Cl 
layer, totally six atomic layers to form this quasi-bulk crystal. To identify different kinds 
of Na atoms, here we mark Nal as Na atoms in the Na-Cl layer, Na2 as Na atoms 
directly touching with the Na-Cl layer, and Na3 as the rest (Figure la). The lattice 
parameters are a=b=3.88 A, c= 33.39 A, a = 90.00°, B = 88.97°, and y = 90.00° 
with atomic positions at Nal, 2e (0.27, 0.75, 0.39), Na2, 2e (0.24, 0.75, 0.54), Na3, 2e 
(0.22, 0.75, 0.69) and Cl, 2e (0.28, 0.75, 0.30). The bond length of Na3-Cl is 2.74 A, 
which is 26% smaller than the Na-Cl bond length in NaCl rock salt phase. It will result 
in an increase of the Coulombic attraction between Na cations and Cl anions. The bond 
length in Na layers, Nal-Nal bond length is 3.74 A, which is slightly (1%) larger than 
the Na bond length of 3.72 A in its metal body center cubic (bcc) phase. The interlayer 
bond length of Nal-Na2 and Na2-Na3 are 3.64 A and 3.96 A, respectively, which are 
slightly (2%) smaller and 6% larger than that in bec phase. 


L3 Éa 


Figure | Structure and Electron Locational Fuction (ELF) of quasi-bulk Na3Cl crystal. 
(a) Side and (b) top view. Na and Cl atoms are presented by blue and green spheres, 
respectively. (c) ELF with an isosurface value of 0.58 e/Bohr*. (d) Two-dimensional 


(2D) slide of ELF on (100) plane. 


The stability of Na3Cl is originated from both the charger-transfer-induced ionic 
interaction and metallic bonding interaction like Na layer in the metallic bec body. In 
the Na-Cl layer, each Cl gains 0.89 e while Nal losses 0.81 e, and charge is closely 
balanced in this layer, like NaCl ionic crystal. Thus, in the Na2 layer, even if Na2 is 
bonded with Cl, there is minor charge transfer of Na2 loses 0.10 e~. In the Na3 layer, 
each Na3 atom gains 0.02 e, similar as the case in metallic bcc Na. If periodic boundary 
condition is applied to this structure, the bulk NaCl; is still stable according to both the 
ionic and metallic interactions, so our result could also be applied in the bulk system. 
Here for the convenience of calculation, we use the six-layer quasi-bulk phase. 

The inner layers of Na show strong metallic properties. The electronic energy band 
diagrams and density of states diagrams of Na3Cl_ P27/m crystals predicted in this work 
are shown in Figure 2a. The contribution of Na to the density of states near the Fermi 
level overwhelms that of Cl. Figure 2b further analyzes the contribution of different Na 
atoms to the density of states near the Fermi level, and the projected density of states 
on Nal and Na2 (inner layers) is much larger than that on Na3 (Na-Cl layer). 

The electronic band structure shows three bands crossing the Fermi level (Ep) in a 


parallel manner contributed by Na. These bands feature a simultaneous occurrence of 
flat bands (T point), corresponding to three sharp van Hove singularities at -2.8 eV, - 


1.9 eV, and -0.7 eV below Er, similar to H3S and LaH0!?”. These electronic structures 
have been shown to be favorable for the formation of stable Cooper pairs, resulting in 


superconductor behavior. 
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Figure 2 Electronic structure of quasi-bulk Na3Cl crystal. (a) The projected band 
structure (left panel) and the projected density of state PDOS (right panel) of 
Na3Cl_P2)/m, the green box represents the enlarged view of the density of states near 
the Fermi energy level. (b)PDOS of ground-state sodium atoms in three different 
environments near the Fermi level. (c) partial charge density (isosurface value of 
0.002e/Bohr3) within leV below Er of Na3C1_P2;/m crystal. 
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Fig. 3. Three-dimensional visualization of energy bands near the Fermi energy level. 
There are three energy bands across the Fermi energy level. Additionally, there is an 
energy band in close proximity to the Fermi energy level. The metallic character of 


Na3Cl_P2//m is attributed to the contribution of the Na atom pair. 


The band structures near the Fermi energy level in Figure 3 indicate the possible 
presence of Fermi surface nesting. The metallic behavior is evident from the bands. We 
can also find nearly parallel bands and contour cut by the Fermi energy level, which is 
consistent with the Fermi surface in Figure 4b and presents the possibility of Fermi 
surface nesting. 

According to the Bardeen-Cooper-Schrieffer (BCS) theory, the calculated 
Electron-Phonon Coupling (EPC) indicates that the superconductivity is due to the soft 
phonon mode. To investigate the superconducting properties of the quasi-bulk Na3Cl 
crystals at normal pressure, the Eliashberg spectral function a@*F(w) and the 
frequency-dependent coupling A(w) were calculated. This result is confirmed by the 


further analysis in Figure.4(a). The coupling constant is calculated by 


> a2F(w) 
A= Zqydgn = 2{ dw ——— 


0 
We can see that the contributions come from the phonon modes near 2-4 THz. It is 
worth noting that the optical vibrational modes corresponding to 2 THz exhibit a 


phonon softening phenomenon, which is reflected by the large projected phonon 


linewidths and a high peak in the Eliashberg spectral function. Phonon vibrational 
modes are involved in the EPC in two different frequency ranges: the low-frequency 
range (< 4 THz), where translational vibrations account for 94% of the total EPC A; the 
mid-frequency branch (>4 THz) is mainly associated with Cl atoms, which contribute 
to 6% of the total à. Figure 4(b) shows the Fermi surface of bulk Na3Cl crystals, with 
a view under the first Brillouin zone on the left and the primitive Brillouin zone view 
on the right. The Fermi surface has three energy bands that cross it, forming a 
multilayered columnar nesting in the first Brillouin zone. The colors on the Fermi 
surface indicate the magnitude of its Fermi velocity, with red indicating a greater Fermi 


velocity. 
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Fig. 4. Calculated phonon properties and Fermi surface of Na3Cl. (a) the phonon 
dispersion curves, phonon density of states (PHDOS), Eliashberg spectral function 
a?F(w) and the accumulated Electron-Phonon Coupling (EPC) constant (œ). (b) The 


Fermi surface of quasi-bulk Na3Cl crystals, with a view under the first Brillouin zone 


on the left and the primitive Brillouin zone view on the right. 


The superconducting transition temperature (TC) is estimated using McMillan- 


Allen-Dynes formula!® 


Wiog —1.04(1 + A) 
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where kg is the Boltzmann constant, u* is the effective screenedCoulomb repulsion 
constant, typically ~ 0.1, à is electron-phonon coupling constant, and Wjog is 
logarithmic average phonon frequency. Phonon softening of acoustic branches is 
observed around the T point. In our quasi-bulk Na3Cl crystals system, A= 0.38 and 
Wiog = 96.63. Using a Coulomb pseudopotential of u* = 0.1, the calculated Te is 0.13 
K. 


Conclusion 

By only adjusting the stoichiometry of Na and Cl as magic ratio Na3Cl, we for the 
first time found the simplest salt crystals daily used can be made an electron-phonon- 
coupled superconductor with a Te of 0.13 K at normal pressure based on first-principles 
calculations. The key to the superconductivity of this quasi-bulk abnormal magic- 
stoichiometric Na3Cl crystal is the cation-rich structure including metallic (Na) atoms 
in the core and hybridization of ionic and metallic bonding, facilitating the electron- 
phonon-coupling for superconductivity. 

It is worth mentioning that the structures and behavior of this quasi-bulk Na3Cl 
crystal could be extended to bulk crystal. Via applying periodic boundary condition 
through the z direction, the stability of this bulk phase is confirmed by the local 
optimization by first-principles calculations. Moreover, the methodology of this bulk 
abnormal crystal constructing may be general to almost all elements, illustrating that 
superconductivity may be possessed by most and even all elements. Notably, its bands 
feature a simultaneous occurrence of flat bands (T point), corresponding to three sharp 
van Hove singularities at -2.8 eV, -1.9 eV, and -0.7 eV below Er, similar to H3S and 
LaHiıo. In the Na3Cl crystals, the relatively large masses of the Na and Cl ions, which 
are not able to produce high frequency phonon modes, result into a very small Te. We 
expect a high value of 7c by using lighter atoms, such as Li or F. Thus, there is hope to 
achieve room temperature superconductivity at normal pressures. 
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